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Abstract 

We use the automorphism group Aut(H), of holes in the lattice Lg = 
A% © A2 © D4, as the starting point in the construction of sphere packings 
in 10 and 12 dimensions. A second lattice, L4 = A2 © A2, enters the con- 
struction because a subgroup of Aut(L±) is isomorphic to Aut(H). The 
lattices Lg and L4, when glued together through this relationship, provide 
an alternative construction of the laminated lattice in twelve dimensions with 
kissing number 648. More interestingly, the action of Aut(H) on L4 defines 
a pair of invariant planes through which dense, non-lattice packings in 10 di- 
mensions can be constructed. The most symmetric of these is aperiodic with 
center density 1 /32. These constructions were prompted by an unexpected 
arrangement of 378 kissing spheres discovered by a search algorithm. 



Introduction 

One of the simplest geometrical constraint problems with strong ties to the design 
of codes is the problem of kissing spheres CO. The object is to pack equal spheres 
— as many as possible — so that each sphere touches a given sphere of the same 
size. The best lower bounds on the maximum number of kissing spheres, in all 
dimensions where this problem has been studied, are derived from integral lattices 
or error correcting codes. Exceptionally symmetric lattices, such as E$ and the 
Leech lattice, account for the relatively few dimensions where the maximum kiss- 
ing number has been established. Conversely, searches for good solutions to the 
kissing spheres problem by an unbiased algorithm tests the scope of the known 
schemes for designing good codes. 
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Figure 1: Results of a numerical search for kissing spheres in 10 dimensions. 
Shown is the distribution of cosines of the angles subtended by pairs of spheres. 
The distribution for 374 spheres has a pair of sharp peaks at cos (0) = (3±\/3)/12. 



This paper comes in response to some unexpected results in the search for 
kissing spheres in 10 dimensions JH. The algorithm used in the search worked 
directly with the two kinds of distance constraints on the sphere centers, Xi : 

ll^ll — ^5 11—^' ^ ~f~ 3 ' 

The only assumption that restricted the search was inversion symmetry (the oper- 
ation X{ — ► — X{ merely interchanges pairs of spheres). A rough characterization 
of a solution is given by the distribution of the cosines of the angles subtended by 
pairs of spheres, 4 cos (6^-) = Xi ■ Xj. Because of inversion symmetry the cosine 
distribution is symmetric about zero. 

For small numbers of spheres there is considerable freedom of movement in the 
solution and the cosine distribution looks quasi-continuous. Figure 1 (left) shows 
a typical output produced by the algorithm for 372 spheres. The number 372 is of 
interest because this is the maximum kissing number of the densest known packing 
in ten dimensions, P\q c fl]. The densest known lattice packing has kissing number 
336. Because there exists a non-lattice packing with kissing number 500 (Piob), 
the quasi-continuous distribution found for 372 spheres is not in itself surprising. 
What came as a surprise was the distribution obtained for 374 spheres. There the 
appearance of sharp peaks, Fig. 1 (right), shows that the sphere arrangement has 
taken a very ordered form. Indeed, with the exception of a very small number of 
spheres, the cosines found by the algorithm are (within numerical precision) all 
contained in the set C = {0, 1/2, (3 db \/3)/12}. The irrational values in this set 
are at odds with the simple rational values generated by integral lattices and error 
correcting codes. 
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A forensic examination of the sphere centers obtained by the algorithm has 
revealed that there is a maximal set of 378 spheres consistent with the cosines C. 
This was accomplished by first finding a suitable coordinate system that reproduced 
the cosines C when the numerical coordinates were approximated by numbers of 
the form a + by/3, a, b £ Q. A basis reduction algorithm was then used to establish 
that the 378 sphere centers form a Z-module of rank 12. 

The main outcome of the forensic examination, and the theme of this paper, 
is a motive for the arrangement of the 378 kissing spheres. In rough outline, we 
first observe that a dense 12 dimensional lattice L\ 2 can be constructed by glu- 
ing together a pair of not exceptionally dense lattices Lg and L4 in eight and four 
dimensions. This gluing is the key to the construction of a dense non-lattice pack- 
ing Qio in ten dimensions. The reduction in dimension occurs entirely within L 4 , 
which possesses a symmetry that decomposes the four dimensional space into a 
pair of invariant planes. Moreover, as this symmetry is non-crystallographic in two 
dimensions, these planes are totally irrational subspaces with respect to L4. The 
packing Q\q results from the projection of L\ 2 into the orthogonal complement of 
one of these irrational planes and explains the rank- 12 Z-module discovered by the 
algorithm. What emerges as a reasonable motive is that the construction of Qiq, 
from the dense L\%, is such that it is able to preserve many of the sphere contacts 
of the lattice packing. 

We will follow the standard terminology of lattices whenever possible [T|. In 
particular, the norm of a vector x is its Euclidean inner product x • x. The automor- 
phism group of a lattice, Aut(A), is the group of unimodular transformations of the 
lattice generators that preserves the Euclidean inner product. Equivalently, Aut(A) 
is the group of real orthogonal transformations on the lattice vectors that permutes 
the elements of A. We use the term glue group for any subgroup of A* /A, where 
A* is the dual of A. The depth of a general point x, in a lattice A, is the minimum 
norm in the set {x — y: y G A, y ^ x}, and is denoted A(x). The number of 
minimum norm elements is written r(x) and equals the kissing number of A when 
x £ A. 



Holes of the lattice A 2 © A 2 © £> 4 

Glue groups of particular interest are those generated by the vertices of the Voronoi 
region, the holes of the lattice. Our construction begins with the lattice 

L s = a 2 e A 2 e £> 4 • 
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Lg has kissing number 6 + 6 + 24 = 36 and for our choice of scale, minimum 
norm 4 and center density 8g = (l/2\/3) 2 (i/8) = 1/96. Each A 2 has two holes 
and the glue group Z3 has a single generator; D4 has three holes and glue group 
Z 2 x Z 2 with two generators. The glue group H generated by the holes of L 8 thus 
has four generators which, for later convenience, we represent as vectors under 
addition modulo 1 : 

= [3 3 3 3] h2 = [I 3 I I] C 1 ) 
h 3 =[± \ 0] /i 4 =[0 £ 

The generator h\ represents a hole in one A 2 component while h 2 represents a hole 
in the other. All four elements, ±/ii and ±/i 2 , thus have depth 4/3 in Lg. Another 
four elements, ±hi ± h 2 , represent holes in both A 2 components and have depth 
8/3. The vectors /14, and h^ + h^ represent the three holes in the D4 component 
of Lg with depth 2. 

The 36 elements of H fall into conjugacy classes with respect to Aut(Lg). The 
quotient of Aut(Ls), with respect to the normal subgroup that acts trivially on H, 
defines Aut(H), the automorphism group of the holes. Aut(H) is generated by 
(i) exchange of the two holes within either A 2 component, (ii) exchange of the two 
^4 2 components, and (iii) any permutation of the three holes of D4. The 2 2 x 2 x 3! 
elements of Aut(H) have an interesting representation in terms of 4 x 4 matrices 
acting on the 4-component glue elements by right multiplication. The matrices 

0-1 " 
-1 
10 
10 _ 

together implement the exchange of hi and h 2 , as well as sending these to their 
inverses (without affecting and h^). Two other matrices generate the 6 permu- 
tations of the three holes of D4. 

0- 10 " 

1- 10 
-1 ' 
1 -1 

Matrix p% cyclically permutes /13, /14 and /13 + /14 (fixing hi and h 2 ), while p 2 
exchanges /13 and /14 (again with no effect on the A 2 holes when combined with an 
appropriate combination of a and P4). 
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Table 1: Properties of the symmetry orbits of glue in L%. The orbit representatives 
h are sums (mod 1) of the generators (Q]); Ag(/t) is the depth and r 8 (/i) the number 
of points in Lg at squared distance Ag(7t). 



For future calculations we tabulate data on the 6 conjugacy classes of the L 8 
glue elements in Table 1. As an example, the entry h = [I | | |] = /ii + /t3 
is the combination of a hole in one A 2 component and a hole in D4. There are 
T$(h) = 3x8 points in Lg at the minimum squared distance Ag(/i) = 4/3 + 2 
from h. The symmetry orbit of h has size 12 and is generated by pzP\- 
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Symmetric glue for the lattice A 2 © A 2 

We now interpret the 36 glue elements, when expressed as 4-component vectors, 
as glue elements of a four dimensional lattice L4. Since we want the glue in L4 to 
include all the automoiphisms the corresponding elements have in Lg, we construct 
L4 so the 4x4 matrices p%, P3, p&, and a are elements of Aut(L±). These matrices 
act on the four generators of L4 by left multiplication. 

We exhibit L4 in terms of the projections of its generators u\, u 2 , U3, U4 in two 
orthogonal planes. These are shown in Figure 2; components in the two planes are 
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Figure 2: Projections of the generators of L4 into two orthogonal planes. 



distinguished by superscripts || and ±. The pair u'^ul^ generates an A2 lattice, as 
does the pair u^,u'l- The same holds for the projection into the other plane. This 
much is consistent with the automorphism generated by p^, which acts as a rotation 
by 120° in the two planes. The pair u 3 ', u\ is rotated by +90° relative to Uj 1 , u 2 \ 
and by —90° in the other plane. If we fix the coordinates so that u\ = uf and 
U2 = v>2, then this implies u 3 ' = — U3 and u\ = —u^. From this we have 

Ui ■ Uj = (uj 1 + uf) ■ (v}j +uj) = ie {1, 2}, j G {3, 4} 

thus identifying L4 as A2 © A2. We choose the scale so that L4 has minimal 
norm 4 and center density 64 = 1/12. The actions of p2 and p4 on the projected 
lattice generators show that these also generate isometries of the two planes: p2 
reflects across a line and p^ rotates by 90°. Altogether, p2, p%, p^ generate the 
two-dimensional non-crystallographic reflection group Go of order 24. The final 
generator of Aut(H), the involution a, reverses the signs of and 114 relative to 
u\ and U2, and is therefore conjugate to an exchange of the two planes. We observe 
that the only point of L4 contained in either plane is the origin. 

The data for the glue elements, as for Lg, again condenses into a table of the 
symmetry orbits. As a guide for the computations in Table 2, we show in Figure 3 
the components of the glue in the A2 plane spanned by u\ and U2 (the diagram for 
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Figure 3: The four types of glue in one A2 component of the lattice L4. 



the other A2 component is identical). There are 12 distinct glue elements in this 
component, and these fall into four families: 



A : 


[0 0] 
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[\ \] 


C : 


[3- a"] 
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D : 


[J i] 


[i I] 


1.3 6J 1.6 3J 1.3 6J L6 6J 



Glue in the family B has depth 1 from 2 points of A2, family C has depth 4/3 
from 3 points, and family D has depth 1/3 from only one point of A2. Taking into 
account both A2 components we obtain the information in Table 2. 

Twelve dimensional lattice packing 

For each element h of the glue group H we have a glue element xs(h) for the 
lattice Lg and a corresponding glue element 2:4 (/i) for the lattice L4. The elements 
xs(h) + X4,(h), h € H, form a glue group for Lg © L4 and, together with Lg © L4, 
form a lattice L12. From the data in Tables 1 and 2 we can check that L\2 has 
minimal norm 4 and compute the kissing number. 

Let x and y be distinct points of L12. The vector x — y belongs to one of 
the 36 conjugacy classes with respect to Lg © L4. A non-zero vector in the class 
[0 0] has minimal norm 4 because it must be a minimal vector of either Lg or 
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h orbit size A 4 (/i) T±(h) 
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Table 2: Properties of the symmetry orbits of glue in L4. A^{h) is the depth and 
T±(h) the number of points in L4 at squared distance A±(h). 



L4; the number of such vectors is 36 + 12, the sum of the kissing numbers of Lg 
and L4. If the conjugacy class h of x — y is non-trivial, then its minimal norm is 
As(h) + Ai{h) and the number of such vectors is T%{h)T±(h). Computations for 
all the conjugacy classes are given in Table 3 and verify that L\i has minimal norm 
4 and kissing number 48 + 96 + 108 + 108 + 288 = 648. The center density of 
this lattice is b\2 = \H\ ds S4 = 1/32. These properties match those of the 12- 
dimensional laminated lattice having the highest kissing number, A^Jf* ID- The 
latter is one of three, equal density lattices produced by laminating one dimension 
at a time, starting with Ai = 2Z and always maximizing the density. L 12 and A^ ax 
were shown to be isomorphic by comparing their Gram matrices 0. 



Ten dimensional non-lattice packings 

Our non-standard construction of = L12 leads rather directly to a non-lattice 
packing that explains the numerical kissing number results in 10 dimensions. The 
construction of the non-lattice packing is closely related to quasicrystal patterns, 
of which the best known example is the Penrose tiling. The focus now shifts to 
the Euclidean space X4 that contains L4. The invariant planes of the the group 
Go provide a natural orthogonal decomposition X4 = X\\ © X±. A general point 
x G L12 has a unique expression of the form x = x\\ + x± + x&, where xs is 
in the space containing Lg. The only role of xg in the following construction is 
through its depth in Lg, which we denote by Ag(x); the depth depends only on the 
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h orbit size A 8 (/i) + A 4 (/i) T%(h)T^(h) number 



\ \ 3 4 32 96 

11114 4 27 108 



l l 

3 3 



4 4 27 108 



5 i 5 l 12 4 24 288 

o 3 o 3 

5 i i 12 6 144 1728 

Table 3: The short vectors of L12 in the non-trivial conjugacy classes of the glue 
group. 



associated glue element. The non-lattice packing is then obtained as follows: 

Construction Q. Define the subset S\q = {x G L i2 : x± G V±}, where V± C X± 
is a bounded domain and specified in detail below. From S10, a packing in 10 
dimensions is given by the set Q10 = x\\ + xg : x G Siq}. 

We will see that if there are kissing spheres at x, y G L12, that is, z = x — y 
is a vector of norm 4, then very often zu ■ zu = z± • z±. Since this implies that 
z\\ + z± + z$ and y/2zu + have the same norm, spheres that are kissing in L\i 
will often have counterparts in Q\q that are kissing as well. Note that if z is in the 
conjugacy class ^\ \ 0],then (y/2 z\\ +z 8 ) ■ {y/2 z\\ +z%) > A 8 (z) = 14/3 > 4, 
and the packing constraint is satisfied for any z\ \ . We derive below a domain V± 
for which we can prove that Q10 is a packing with minimal distance 2. 

As a preliminary step we identify and characterize the minimal vectors of Lu 
that arise in construction Q: 

Lemma. A non-zero x G L12 that satisfies the inequalities 

A 8 (x) + 2x\\ -x\\ < 4 (2) 
x± ■ x± < 8/3 (3) 

is a minimal vector of L 12. 

Proof. Combining (2) and (3) we obtain a; ■ x < 14/3 + ^ Ag(x). Table 1, together 
with A 8 (x) < 4 as implied by (2), limits the possible values of A 8 in the bound 
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representative X H ' x || x± ■ x± 



1111 4 4 

3 3 3 3 3 3 

11 _2 1 4 2 4.2 

3 3 3 3 3 ^3^^ 

11 1 _2 4|_2_4 2_ 

33 3 3 3 ^3 3 ^ 

Table 4: Data on the three orbits with respect to the reflection group Go of the 4 x 9 
minimal vectors of L^i with A4 = 8/3. Orbit representatives are specified in the 
same 4-component basis used to define the glue group. 



on x • x. Recalling that the possible norms of the integral lattice L\i = A^ ax are 
4, 6, . . ., we see that only the case A§(x) = 10/3 is unresolved. In this case (2) 
and (3) imply xm • < 1/3 and X4 • X4 < 3, where X4 = xu + x±. Without 
loss of generality it is sufficient to examine just one element of the symmetry orbit 
of the glue, say h = [g g g g] • There are just five solutions to the inequality 
X4 • X4 = || Yli=i { n i + g) u i\\ 2 < ^ for integers n; L . These are where all integers 
are zero, for which X4 • X4 = 2/3 and x is minimal, or only one of the integers is 
nonzero and has the value —1. For the latter one obtains x\\ ■ x\\ = (4 ± \/3)/3, 
both of which are inconsistent with X|| ■ x\\ < 1/3. □ 



The domain V± is constrained by the symmetry orbits of the L12 minimal vec- 
tors, specifically their projections into X±. Since Go acts as the reflection group 
of order 24 on the set of projected minimal vectors, the size of the orbits is always 
a multiple of 12 (since x± is non-zero for a minimal vector). Consulting Table 2 
we see, for example, that there are 3 x 4 minimal vectors with X4 • X4 = 2. This 
single orbit with respect to Go is also an orbit of Aut(H) and is therefore fixed by 
the involution a that exchanges Xu and X±. We conclude that for minimal vectors 
x in this conjugacy class we have xm ■ xm = x± ■ x±. The equality of projected 
minimal vector norms in X\ | and X± holds for all the classes except the class with 
Ag = 14/3, where this property will not be relevant, and the class with A§ = 4/3. 
The latter contains 4x9 vectors which form three orbits with respect to Go; data 
for these are given in Table 4. 

Only the second orbit in Table 4 requires attention, as it violates the property 
llv^xull > 1 1 xi I + xj_||. Avoiding such minimal vectors is one of the roles of the 
domain V± . More specifically, we construct V± such that if y± , z± £ V±, then 
xj_ = y± — z± is never one of the 12 forbidden vectors in the orbit of (uj- + 
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Figure 4: The domain Vj_, a regular dodecagon, used in construction Q. Trans- 
lations associated with pairs of opposite edges (arrows) are the forbidden vectors 
of the construction. V± includes half the boundary (dark), and in particular, six 
dodecagon vertices. 

U2 — 2 u% + )/3. Up to translation in X±, there is a unique dodecagon with the 
property that its opposite edges are translated by vectors in this set. It is possible to 
include half the boundary of the dodecagon in the definition of V± , and we make 
the definite choice shown in Figure 4 in what follows. V± has diameter y/8/3, 
area |Vj_| = 2, and happens to coincide with the shadow of the Voronoi domain of 
L4 = A2 © A2, although the significance of this is not clear. 

Theorem. Construction Q, with V± as specified above, produces a sphere packing 
Qiq with minimum distance 2. 

Proof. Consider x, y € S\q and let z = x — y. We wish to verify that \\y/2 zu + 
zs 1 1 > 2. From the construction of S±o we know that z± • z± < 8/3 (the diameter 
of Vj_). The theorem is proved if we can show that ||-v/2^|| + ^8 1 1 < 2 leads to a 
contradiction. But with this statement both hypotheses of the lemma are satisfied, 
thereby establishing that z is a minimal vector. From our exhaustive analysis of the 
minimal vectors we know (again omitting the irrelevant class with A§ = 14/3 > 4) 
that \\y/2 z\\ || > ||z|| + z±\\ with one exception: the minimal vectors whose projec- 
tions into X± coincide with one of the forbidden vectors. But V± was designed so 
this is impossible, hence 1 1 y/2 z\ \ \ \ > \ \ z\ \ + z± 1 1 . Moreover, since z is minimal, this 
shows ||\/2;z|| + z%\\ > \\z\\ + z± + z&\\ = 2, a contradiction. □ 

We know that the packing Qiq is aperiodic because its projection into X\\ is 
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Figure 5: Sphere centers (dots) of the packing Q\q projected into the plane Xu form 
an aperiodic pattern, a region of which is shown here. Pairs of kissing spheres in 
the 10 dimensional packing may have any of the five projections shown. 



a quasicrystal. A representative region of the projection is shown in Figure 5, 
where each point represents an Lg lattice translated by one of 36 glue elements. 
To set the scale, the projections of various pairs of unit-radius spheres are also 
shown. Even though the sphere projections overlap, glue translations of the Lg 
lattices guarantee that the spheres form a packing. There is a smallest separation of 
projected sphere centers, and spheres in this relationship amply satisfy the packing 
constraint because the associated conjugacy class of the gluing has A§ = 14/3. 
Many spheres are actually kissing; this is the case for those depicted in the figure. 

The center density of Qio can be computed using a formula developed in the 
study of quasicrystals |U[5). Restricting to just one of the conjugacy classes of 
glue, and omitting the dilation by the density of centers in Xu is equal to the 
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Figure 6: The projected sphere centers in form a quasicrystal tiling comprising 
triangles, squares, and 30° rhombi. 



product of the center density of L4 and the area of the domain V± whose shadow 
selects the subset of L4 used in Qiq: Su = \V±\ 64. Including all the conjugacy 
classes multiplies this by \H\ =36 and the dilation by \/2 diminishes the density 
by 2. To get the density in 10 dimensions this density in is multiplied by the 
center density of L§: 

Because Qiq is a non-lattice packing, the sphere centers fall into distinct classes. 
This is evident in the quasicrystal pattern of projected centers in X\ | , where there 
is a special class of centers that is never in the shortest-separation relationship with 
other centers. When centers in this special class are connected by edges, the re- 
sult, shown in Figure 6, is a tiling of triangles, squares and 30° rhombi. Alternate 
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tilings formed by these three tiles, utilizing the same gluing scheme for the Lg lat- 
tices, also correspond to valid packings (details omitted). The densest of these is 
obtained by the tiling that only uses the triangle, the densest tile of the three. The 
improvement of the density, to the value 7/(96 + 64\/3) ~ 0.0338, comes at the 
expense of symmetry. Whereas the quasicrystal pattern has point symmetry group 
Go of order 24, the point group of the triangular tiling is only of order 12. The 
densest known packing in 10 dimensions has center density 5/128 « 0.03906 OX 

We conclude by returning to the results of the numerical experiment that prompted 
this investigation. The arrangement of 378 kissing spheres discovered by the search 
algorithm [2] coincides with the arrangement obtained from the packing Qiq when 
the dodecagonal domain V± is given a singular centering (translation) in X±. At a 
singular centering one Lg lattice projects to the exact center of V±, thereby making 
the six vertices of V± available to the packing (see Fig. 4). Relative to the do- 
decagon center, the six vertices have glue in the orbit of [| | 0] with Ag = 8/3 
and A4 = 4/3. From Table 1 we see that rg = 9 spheres, in each of the Lg lattices 
that project to these vertices of V±, make contact with the central sphere. There is 
also a set of 12 Lg lattices that project to a regular dodecagon within V±. These 
are in the orbit of [§ § | §] with A 8 = 10/3, A 4 = 2/3 and r 8 = 24. The X\\ 
counterpart of this relationship are the many examples of complete dodecagons en- 
circling tile vertices in Fig. 6. Finally, the Lg lattice that projects to the dodecagon 
center makes 36 contacts with the central sphere (the kissing number of Lg). The 
net kissing number is thus 6x9 + 12x 24 + 36 = 378. 
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